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Abstract 



o 
o 

(N 

^ ■ Nonperturbative self-energy of a bound quark is computed gauge-invariantly in 

the framework of background perturbation theory. The resulting Am^ is negative 

in 



and is a universal function of string tension and current quark mass. The shift of 

2 



psj ■ hadron mass squared, M , due to Am^ is negative, large for light quarks, and solves 

<*' , the long-standing problem of Regge intercepts in relativistic quark models. 

^ ■ 
^ ■ 

in : 1 Introduction 
o 

The quark self-energy (QSE) operator for a quark propagating through the vacuum back- 
ground field can be introduced in a gauge-invariant way, when quark is a part of a white 
Ph! object. Then one associates QSE with a part of the total energy of the bound state which 

JlLii does not depend on characteristics of other constituents, but depends on the given quark 

^ ' motion. In contrast to the perturbative QSE which is gauge and renormalization-scheme 

dependent [||] one can introduce the nonperturbative QSE in the framework of the Back- 
. ground Perturbation Theory (BPTh) @, |[, in which case QSE can be expressed through 

! the gauge invariant field correlators (for a general formalism of field correlators see |0, P ) • 

I It is useful to separate the nonperturbative contributions to QSE into two parts; one, 

coming from the color quark charge interaction with the nonperturbative background, 
and another created by the color magnetic moment of the quark (analogs of diamagnetic 
and paramagnetic contributions). 

It was shown in |5[ that the quark charge part is responsible for creation of confining 
interaction hence is nonlocal and should be attributed to the qq interacting kernel. Only 
small distance-independent corrections can be separated out of this part and associated 
either with QSE or with qq interaction. 

The color-magnetic moment part of QSE is instead large and independent of the string 
and antiquark motion (in the limit of small correlation length Tg of field correlators). It 
gives the dominant part of QSE and has specific negative sign, decreasing the hadron 
mass. In what follows we shall compute only this part and shortly discuss corrections 
from the charge part in the conclusions. 

The next step is to calculate the contribution of QSE to the hadron mass and Regge 
trajectories. To this end one can use the Feynman-Schwinger (world-line) representation 



and following from it the Hamiltonian einbein technic, introduced in |1lO| , [TT[] (for 



reviews see lectures and recent paper [jT3[)- As a result one obtains a negative 
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self-energy contribution from each quark and antiquark depending only on string tension 
and current quark mass. 

Taken as a correction to the total hadron mass, this selfenergy contribution is favoured 
from phenomenological point of view. First of all it explains sign and magnitude of 
negative constant which is always added to the relativistic quark Hamiltonian in most 
existing models. In a similar way is resolved the old problem of Regge intercepts J(M^ = 
0), when intercepts come out too low when computed through string tension without 
introduction of that phenomenological constant. Some examples for meson and baryon 
masses are given and the net effect for the Regge slope is discussed. The nontrivial 
character of the hadron mass shift obtained below is that it corrects the Regge intercept 
and keeps at the same time the value of the Regge slope at the standard string value I'kg. 
(One should note that simply adding a constant to the hadron mass, as it is usually done 
phenomenologically, would spoil the universality of Regge slope). 

The plan of the paper is as follows. In section 2 the quark Green's function in the 
framework of BPTh is introduced and the world-line integral is used to separate contri- 
butions of color charge and color magnetic moment of the quark. In section 3 the lowest 
order (Gaussian) contribution to QSE is written down and expressed through the Gaus- 
sian correlators. In section 4 self-energy corrections to the Hamiltonian are derived and 
estimated, and the problem of Regge intercepts is discussed, while section 5 is devoted to 
conclusions and outlook. 



2 Quark Green's function in vacuum background 
fields 

Gluons play two different roles in QCD, namely, they create gluon condensate and string 
between quark and antiquark on one hand, and also propagate in the nonperturbative 
vacuum as valence gluons. Therefore it is convenient following p|, || to separate the total 
gluonic field into background NP field and field of valence gluons a^, 

= 5^ + (1) 

and using the 'tHooft identity to write the QCD partition function as 

with N' = N J DBri{B). Here 1]{B) is an arbitrary measure and can be put 1]{B) = 1 or 
chosen in such a way as to improve convergence of perturbation series in ga^. For what 
follows the value of r]{B) and explicit details of averaging over fields -B^ are not important. 
According to (^ one can write any vacuum average as 

{F{A))A = mB + a)))B,a (3) 

The Euclidean quark Green's function Gq{x,y) can be written using the Feynman- 
Schwinger (world-line) Representation (FSR) 10-0 as 

G,(x, y) = (m, + D)-J = (m, - D)^{ml - D^)-^ = 
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roo rx rs 

= {mg-D)^ ds{Dz)^ye~^PAexp{ig A^dz^)PFexp{ gaFdr), (4) 

Jo Jy Jo 

where Pa, Pf are ordering operators of matrices and F^,^ respectively, aF = a^yF^^ = 

( ctB ) ' ^^^^ ^f"" ^ Tii^t'^'' ~ ^^^z^)' K = ^ls + l /o z^dT. 

In (^) the quark moving along its trajectory z^{t) is interacting with external total 
field Afj^ by its color charge (the first exponent) and by its color magnetic moment (the 
second exponent in (^). As we shall see it is the second interaction which yields the 
most important effect in creating the nonperturbative mass shift of the quark, and we 
shall calculate it in terms of field correlators. One should also note that the term ga^ in 
D = J^iidfj, — iga^ — igB^) produces usual perturbative quark mass correction (including 
anomalous dimension for the mass evolution) and was investigated thoroughly especially 



for heavy quarks , (see and references therein), and we shall not enter in discus- 
sion of this topic. Instead we shall concentrate on the contribution of nonperturbative 
background B^. 

The Green's function (^ is not gauge-invariant. To define quark mass corrections in a 
gauge-invariant way, one must use the Green's function of a white object - of a meson or 
baryon. In the former case one can write [|15] (in the fiavour-nonsinglet case and neglecting 
sea-quark determinant) 

GM{x,y) = {trTiGq{x,y)T2Gg{x,y))B = 

^ /"CO roo 

= (tr{ri(mq, - D) dsiDz ds2Dze-^'-^^<l>„{x,y)x 
Jo JO 

xT2{mg^- D)<^a{y,x)})B. (5) 

Here $(j(x, y) is the product of the last two exponents on the r.h.s. of (|^), where the 
field Afj^ is replaced by 5^ (setting = 0): $o-(x, y) = PbPf exp[i5f B^dz^ + Jq gaFdr]. 

One can note that the integrand of (^ contains the closed Wilson loop with insertions 
of operators D (twice) and operators aF (infinitely many times). Therefore the whole 
construction in (^ is gauge invariant. The sign of trace in (^) implies summing over 
Lorentz and color indices, and rj(z = 1,2) stands for current vertices, Fj = 1,75,7^, ... 

3 Quark self-energy correction due to nonperturba- 
tive background 

Consider the FSR for the quark Green's function (^). In the 2-nd order of perturbative 
expansion it can be written as 

^ roo rs rTi 

Gq{x,y) = {mq — D) J ds j dri j dT2e~^ {Dz)xud'^u{Dz)uvd'^v{Dz)yyX 

{igA^{u)u^ + ga^yF^y{u)){igAy{v)vy + g^xaFxaiv)), (6) 
where we have used the identities 

{Dz)^y = {Dz)^ud'^u{Dz)uvd'^v{Dz)^y, (7) 
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roo rs rTi roo roo roo 

/ ds dri / dT2f{s,Ti,T2)= / ds dri / rfr2/(s + n + rg, n + ra, Ts) (8) 
Jo Jo Jo Jo Jo Jo 

At this point one can expand the quark Green's function only in color magnetic moment 

interaction (cxF), which is useful when spin-dependent interaction can be treated pertur- 



batively, as it is in most cases for mesons and baryons ( the exclusions are Goldstone 
bosons and nucleons, where spin interaction is very important and interconnected with 
chiral dynamics). In this case to the second order in {(tF) one obtains 

^ /"oo foo roo ^ 

Gf\x,y) = {m„ - D) ds dri dT2e-'^'^^'+^^+^^^-^''-^^-^^{Dz)^^x 
^ Jo Jo Jo 

X <^{x,u)g{aF{u))d\{Dz)uvHu,v)g{aF{v))d%{Dz)^y. (9) 

In another way expanding Gq = {nig — D){m'^ — — gaF)^^ in powers of {gcF) it 
can be written as 

Gf^x^y) = K - D){ml - Dl&u g{aF{u)){ml - D^-^Sx 

^g{aF{v)){ml-Dl)-^. (10) 

One can easily find the correspondence between (|]) and ( |T0|) since (m^ — -D^)^^ is the 
scalar Green's function having the FSR as follows 

/•oo 

y) ^ {ml - Dl)-l = ds{DzUe-''<^,{x, y). (11) 

Now keeping terms g{(yF) up to the second order one can find from ( [T0| ) and (^ 
(leaving aside for simplicity the factors {rrtq — D) in (^) that one has in (^ an expansion 

tr{W)B + tr{PFWgaF{u)gaF{v))B + ... (12) 

where operators gaF are inserted inside Wilson loop in the properly ordered form due to 
the operator Pp, keeping gauge invariance. 

Using cluster expansion for W ^ the last term in (|T2|) was computed in the Appendix 



of with the result 

{trPFF^,{u)Fpx{v)W{C)) = tr{[{F^,{u,xo)Fpx{v,xo))- 



-g J daai3{w){F^^{u,Xo)Fc,f3{w,Xo)) J da^s{z){F^x{v, Xo)j5F{z, xo))]{W{C))} (13) 

Here F{x,xo) = $(xo, a;)F(x)$(x, xq), with ^{x,y) = Pexpig jy Bfj,{z)dz. It is conve- 
nient in (|13D to choose the point xq (note that {W{C)) does not depend on xq) on the 
line between points u and v, since in this case the first correlator in ( P^ takes the form, 
studied in [g 

g^{F^,{u)^u,v)F,x{vMv,u)) = 

= ^ab{{S^,pSux - Sf,x - SpxSup)D{u -v) + ^[d^{hp6^x - hx6up) + 

+ d^{hx5pp- hpSpx)]Di{u-v)}, hf, = u^-Vp, (14) 
where lab is the unit operator in color space. 
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In what follows we shall keep only the first term in Eq. (|I^) in the form ( (|T^) , since the 
second term describes interaction of the quark spin with the string world sheet between 
the quark and antiquark and cannot be clearly associated with self-energy (actually this 
term contributes to the spin-orbit interaction between quark and antiquark, see discussion 

0). 



m 



Multiplying (|T^), (0) with spin operators one has 

g\<yF{u)^{u, v)aF{v)^{v, u)) = iabia^{8D{u - v) + ^d,{hMu - v))}. (15) 
Here lap is the unit operator in bispinor indices. 

iFiom ([T0| ) one can see that the following combination occurs playing the role of quark 
selfenergy correction 

A = J d\u - v)g^aF{u)^{u, v)aF{v)^{v, u))G{u, v) (16) 

where G{u, v) = ("^g — -0^)^,1, is given in (11). Since G{u, v) depends on the interaction of 
quark with the string world sheet (through = d^ — igB^), then in general the correction 
(|I6|) cannot be unambiguously attributed to QSE. However, at this point one can use 
the important fact, that D{u) and Di{u) are short-ranged functions with the gluonic 



correlation length Tg 0.2 ^ 0.3 fm from lattice calculations in the cooled vacuum |T7 
and even smaller, if one extracts Tg from the gluelump masses |]l8l . Therefore integration 
in ( p!6D is limited to small \u—v\ and in this region one can replace G{u, v) by Gq{u, v) (since 
Go yields singular terms ~ ( ^ \2 i while background field corrections give nonsingular 

[U V) 

contributions of the order of Tg (F^) . Thus one can expect the accuracy of this replacement 
to be of the order of g'^Tg{Ef) ~ 10%). Replacing G by Go in ( |16D and using for D, Di 
exponential form found on the lattice [0, one obtains in this approximation from (|I6|) 

Ao = 1 d^w{8D{0)e-'\^^ + ^D^{0)^,{w,e-'\^\)}x 

S^^- (17) 
It is convenient to define an integral {Ki{x) is the Mc Donald function) 

/•oo 

(p{mg,6)= w'^dwe-^'"Ki{mgw) (18) 
Jo 



mqKi{mq 


w\) 


47r2 


w 





which has analytic form for 5 > m, 



q 



(5^ — m^)^/^ TUg rng{6' 

In the opposite case, S < mg, one has 



ip{mq,d) = -— — — — ^arctan + — — — . 20 
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Two limiting cases are of interest for what follows 



^{mq,5)[m„^Q= (21) 



^(m,. = ^3-1^ + 0(^1). (22) 



In a similar way one can define the integral in front of -Di(O) in Eq. (|T^, 

d_ 

86' 



ipiinig, 5) = (4 + 6—)(p{mg, 6) 



In _ .^^ ; > 



(5^ — m'^y/'^ niq mq{6'^ — m- 



15"^^ a..ta.V^+^"'t^"'n:f^(^<>n.), (23) 



The limiting values of are 

2 

<^l(?^g,5)|m,^0 = — n (24) 

^i(m„5)U,_ = --^ + 0(^-j. (25) 

Comparing the expansion of (m^ + Am^ — D^)^^ with (p!OD one can define the QSE 
for niq 

Ami = = -^(41^(0) + ^1^1 (0)). (26) 

For purely exponential correlator D{x) one can connect -D(O) with string tension a; in the 
Gaussian approximation, which was checked recently using Casimir scaling arguments to 
be accurate within few percents [|l^, one has 

As a result one has for the self-energy if rriq —>■ 

Lattice measurements 0] give ^ < ^, and later we shall omit this term. We note at this 
point that Eq. ( pHf ) is universal in the sense that it does not depend on the exponential or 
any other form of correlator in (^), since the integral (0) for niq —>■ assumes the form 
of integral (P7| ) for string tension. 

For nonzero values of rriq one has instead of (|28| ) 

Amq{mq} = - —V, with r] = . (29) 



TT 



For example, for the strange quark with rris = 0.175 GeV one has 1] = 0.88, while for 
the c-quark, rric = 1.7 GeV, rj = 0.234 and finally for the 6-quark, with rrib ~ 5 GeV one 
obtains r] = 0.052. 
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4 Self-energy corrections to Hamiltonian and hadron 

masses 

Using expression (H) for the meson Green's function and exploiting the einbein formahsm 



i=l 



+ 4!/'4§,W'4^,;3 (30) 



and 



10| , pA| one can construct the relativistic Hamihonian. We introduce the squared quark 
masses with corrections due to QSE obtained in the previous section, — > + Am^ 
and arrive at the expression 

H = Ho + AHq + AHs, 

° ^1 2^ 2/2 2[^i(l-C)2 + /i2C' + /o^W-C)M/5)] 

W) 

AH, = E — (31) 
Here we have defined the reduced dynamical mass p, and the parameter (: 

and and mg{2) are current quark and antiquark masses, renormalized at the scale 

around 1 GeV, while fii, fi2 and are einbein functions which should be found from the 
minimum of this Hamiltonian. For the discussion and practical use of this Hamiltonian 
including spin-depending part AHs, beyond the original papers |jll| see also lectures 
and recent paper 



It is important that the dynamical masses Wi,yU2 are found from the stationary point 
analysis of the Hamiltonian and these values 1^1,1^2'^ can be used to estimate self- energy 
corrections to the Hamiltonian and its eigenvalues. 

Hence the eigenvalues of H will be shifted by an amount read off from ( ^l]) 



Ai^™ = --f4T + 4wV (33) 



i.e. from each quark the shift is equal approximately to jgy. Taking into account that 

TTfJ.\ 

the lowest eigenvalue M^^^ is equal (without color Coulomb interaction) 4/x*^°) |T^, |11|, |T3[, 
one has for equal masses fif^ = /i2°'' = /i^^-* 

M(o) {L = nr = 0)+AM = 4/i(o) - (34) 



one can 



In a similar way one can write the shift of the baryon mass. Using [^T], |22 
write 

2cr 3 ]^ 

AHq{baryon) = ^ (35) 

^ 1=1 K ^ 
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where are to be found again from the minimum of baryon Hamiltonian In 
the lowest approximation (color Coulomb and spin interaction neglected) and for light 
quarks one has /if^ = and M*^°) = 6fi^^\ so that the resulting corrected mass is 

M = M(° + AH„ = 6/i(°) - (36) 



The same type of correction occurs for hybrids, since AHg has the same form (33) 



as for mesons, however the values of nf*^ and are to be found for the total hybrid 



Hamiltonian (see |]12|) containing in addition dynamical mass /i^ of the valence gluon. 

We now come to the one of the most important problems of hadron spectrum, which 
to the knowledge of the author had before no solutions - the problem of Regge intercepts. 



It was shown in |2^, |26|] that the Hamiltonian Hq ensures the string Regge slope and 
the square of its mass eigenvalue is M'^{L) = 2'KaL + Cq, where Cq is positive and too 
large (for = 0, Cq ~ llcr = (1.4 GeV)^ when Coulomb and hyperfine interaction is 
not taken into account), hence the Regge intercept 

J(M = 0) = L + 1 = 1- -^ (37) 

27ro" 

is negative (-0.73), while for the experimental p trajectory it is positive and near +0.5. 
Therefore one needs such mass corrections which would shift Regge trajectories upwards 
in the Chew-Frautschi plot (J(M^) vs M^) not changing its slope. 

As we shall see now, the correction Aifg (|30D has exactly this property. To start we 



take the corrected mass ( p6| ) and consider its square 



M2 = (m(°) + AMf = (4/i(°) - = {M^^^f - ^ + 0(AM)2) (38) 

Thus to our accuracy (first order in AM) mass squared is shifted by a constant not 
depending on quantum numbers, — and if (M*^"))^ corresponded to a linear Regge 
trajectory, the same will be true for the shifted trajectory with the same slope. 

If one calculates M '^''^ = from the simplified Hamiltonian, obtained from (|3lD by 
putting L = and pi — >■ p^(this is very close to the relativistic quark model Hamiltonian 
Hrqm = 2i/p^ + + ar) then the slope of M*^°) would be 80" instead of expected string- 



like 27rcr [jTO[, and this slope is kept intact by the shift (|38D . To resolve this inconsistency 
one should take into account the term with in (^Tj) which contains the moment of 
intertia of the rotating string and changes the slope to the correct value 2na (pT]). 



Exact calculations with the Hamiltonian (|3l|) require numerical solutions yielding the 



linear Regge trajectories with the slope very close to 2'7ia (see [p5| , |26| for results and 
discussions). Below we shall consider the effect of string rotation as a correction 
[]1T[ and write for the whole mass squared 

M2^(4/iW + =^ + AM(i))2 (39) 
with Ami = - V and AM^^) = -f "'^/(^ffl . Denoting 

(4/i(°) + -j^y = 8aL + Co, (40) 
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one obtains to the lowest order 



The Regge slope can be computed from (^) both for L = 1 and for L — » oo, 
which yields ^ = f a(L ^ oo), = 6.5a(L = 1). 

One can see that both slopes are rather close to the string one, 27ra, while the self- 
energy term Am^ contributes only to the constant Cq, shifting it towards the experimental 
value of m^. Thus QSE indeed can solve the problem of Regge intercepts to a reasonable 
accuracy. To get an idea of the magnitude of AHg, consider the p - meson, corresponding 
to the solution of the Hamiltonian (|30|) with L = n,. = 0. Taking cr = 0.18 GeV^ one 
obtains 0, |12| /xf^ = = /i^o) = 0.352 GeV and M^o) = 4^(o) _ ^ q^^^ ^^lile 
(|3|) yields AHg{meson) = -0.65 GeV, with the total mass M = M(°) + AHg = 0.76 
GeV. This last figure is close to the experimental p-meson mass, but our result should be 
corrected by hyperfine spin term and Coulomb interaction which tend mass in different 
directions with a net shift around or less than 0.1 GeV. 

A similar situation occurs for baryons. Taking a = 0.15 GeV^ (see and for the 



phenomenologically motivated choice and its theoretical justification), one obtains from 
m for light quarks /i^") = 0.957^0", and = Qp^^l As a result one obtains for A 



mass (without spin and Coulomb correction) 

mJj"^ + AM = (2.22 - 0.77)GeV = 1.44GeV. 

For strange baryon, Q, one has instead for mq{s) = 0.175 GeV the dynamical mass 
Us =0.415 GeV O and finally 

M^°\n) + AHg = (2.46 - 0.607)GeV = 1.85GeV. 

These examples are only to illustrate the order of magnitude of resulting corrections. 
We leave the problem of realistic meson and baryon mass calculations with spin and 
Coulomb interaction included to future publications. 



5 Discussion and conclusions 

In most model calculations of hadron masses done heretofore the absolute values of masses 
have been defined up to a constant, varying from one family of hadrons to another. Also 



in our previous calculations |T^, |T3[ this constant was introduced and associated with the 
quark self-energies, however the latter were defined phenomenologically [0 to be around 
-0.25 GeV per light quark, but not derived from theory. 

In the present paper we have looked into problem of quark self-energies and absolute 
values of hadron masses, and calculated both in terms of the only scale of nonperturbative 
QCD - the string tension a. 

As a result we have succeeded in several respects. Firstly, all our theory is based on 
quark current masses, assumed to be renormalized at the scale of around 1 GeV, where 
NP effects do not yet practically enter. These masses should be associated with the pole 
masses, introduced in case of heavy quarks (see e.g. |I^]). 
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Secondly, QSE corrections yield a numerically reasonable shift of masses downwards, 
and for light quark this shift is indeed around -0.25 to -0.3 GeV, as was required by 
phenomenology |1T3| . 

Thirdly, and this is most important featurefrom theoretical point of view, this QSE 
shift does not spoil the correct Regge slope and at the same time resolves the old problem 
of Regge intercept which was too low in all previous calculations (if negative constant is 
not introduced). 

Let us now discuss the nature of this correction and compare it to another possible 
approach. First of all, the negative sign of QSE is evidently connected to the paramagnetic 
mechanism, which is at work here. One should note that the quark charge interaction, 
present in (H) in the form of gAf^, produces effect of another (positive) sign to the quark 
mass squared, and the net result is bounded from below, since in (m^ — £)^)~^ the operator 
—D^ = —D^ — g{aF) is nonnegatively defined, and vanishes only for quark zero modes. 
However in the confining phase the charge part gAf^^ produces for the white system of [qq) 
the confining interaction ar and belongs to the whole worldsheet of the string and cannot 
be associated with QSE. Instead the part g{(TF) contains the local correction (in the limit 
of small Tg) g'^{{aF){aF)) , which can be treated as QSE and was computed in the present 
paper. One important conclusion is that actually QSE should be considered finally as a 
correction to the total hadron mass, and not to the mass of the quark, since the latter 
has no definite meaning inside the hadron. Therefore one should not give much sense to 
the fact that for light quarks the total quark mass + Am^ is negative -one should 
have in mind that this contribution is only a part of the Hamiltonian (^) and should be 
computed as a correction to its eigenvalues. (Note that one cannot include AHg in Hq to 
find jjf'^ from stationary point analysis, since this procedure would go beyond the lowest 
order in AHq, while AHg was itself computed to the lowest order only). 

To get more understanding of negative QSE correction to the quadratic (in Dirac 
operator D) Hamiltonian, one can consider masses of heavy-light mesons Mhl = Mq + e 
computed in two different approaches, namely to compare eigenvalues of this quadratic 
Hamiltonian in the case when one mass, m^^^ = Mq, is very heavy, to the eigenvalues 
of the corresponding Dirac equation with linear confining interaction. Taking the latter 
from [p7l B3] one has for lowest Dirac eigenvalues. 



ei!'^ = 1.62v^ (42) 
while for the Hamiltonian Ho(m^^^ oo) one obtains ||13| 

ei"'^ = 2.26v^. (43) 

One can see the difference of about 0.3 GeV, which is very close to the mass shift due to 
QSE. This is not surprising since in Eq^^ the corresponding QSE correction is absent, while 

for this correction is present and is given by -g^- Therefore working in the linear 
Dirac formalism like that developed in p9| one should not include the aforementioned QSE 



correction and the Regge-intercept problem does not appear. Finally it is interesting 
to note, that self-energy correction similar to the one considered in the present paper 
appears in the 1-1-1 QCD, where in 't Hooft equation each quark obtains a negative 
contribution to the mass squared equal to — 2o"/7r, i.e. one half of self-energy correction 
computed above in 3-f 1 QCD. 
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Summarizing, we have calculated the nonperturbative QSE correction explicitly and 
found the shift in hadron masses and Regge intercept which is strongly favored phe- 
nomenologically. 

Results of more detailed analysis for both light and heavy hadrons are to be published 
elsewhere. 

The author is greatly indebted to A.M.Badalian and A.B.Kaidalov for many stimu- 
lating discussions. This work was financially supported by the RFFI grants 00-02-17836 
and 0015-96786. 
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